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Abstract. We report the analysis of an XMM-Newton observation of the close binary HD 159176 (07 V + 07 V). The observed 
Lx/Lboi ratio reveals an X-ray luminosity exceeding by a factor ~7 the expected value for X-ray emission from single O-stars, 
therefore suggesting a wind-wind interaction scenario. EPIC and RGS spectra are fitted consistently with a two temperature 
mekal optically thin thermal plasma model, with temperatures ranging from ~2 to 6 10 6 K. At first sight, these rather low 
temperatures are consistent with the expectations for a close binary system where the winds collide well before reaching their 
terminal velocities. We also investigate the variability of the X-ray light curve of HD 159176 on various short time scales. 
No significant variability is found and we conclude that if hydrodynamical instabilities exist in the wind interaction region of 
HD 159176, they are not sufficient to produce an observable signature in the X-ray emission. Hydrodynamic simulations using 
wind parameters from the literature reveal some puzzling discrepancies. The most striking one concerns the predicted X-ray 
luminosity which is one or more orders of magnitude larger than the observed one. A significant reduction of the mass loss 
rate of the components compared to the values quoted in the literature alleviates the discrepancy but is not sufficient to fully 
account for the observed luminosity. Because hydrodynamical models are best for the adiabatic case whereas the colliding 
winds in HD 159176 are most likely highly radiative, a totally new approach has been envisaged, using a geometrical steady- 
state colliding wind model suitable for the case of radiative winds. This model successfully reproduces the spectral shape of the 
EPIC spectrum, but further developments are still needed to alleviate the disagreement between theoretical and observed X-ray 
luminosities. 
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1. Introduction 

HD 159176 is a relatively bright (m v = 5.7) double-lined spec- 
troscopic binary (Trumpler 1930) in the young open cluster 
NGC 6383. The system has been well studied in the visible and 
UV wavelengths. Conti et al. i 19751 1 derived an orbital solu- 
tion that was subsequently improved by Seggewiss & de Groot 
dl976> . Lloyd Evans i 19791 ) and most recently by Stickland et 
al. (Q993). The binary has an orbital period of 3.367 days and 
consists of two nearly identical O-stars in a circular orbit. 

Conti et al. suggested that both stars are 07 stars that 
have evolved off the main-sequence and nearly fill up their 
Roche lobes. However, Stickland et al. argued that the stars 
were probably not evolved. Although the system does not 
display photometric eclipses, Thomas & Pachoulakis (1994) 
reported ellipsoidal variability with an amplitude of about 
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0.05 mag in the optical and UV wavebands. These light curves 
were analysed by Pachoulakis (1996) who inferred radii of 
both stars of order 0.25 a, where a is the orbital separation. 
According to these results, the stars do not fill up their critical 
volume and are thus not very deformed. 

Several observations point towards the existence of a wind 
interaction process in HD 159176. For instance, the optical 
spectrum of the system displays the so-called Struve-Sahade 
effect, i.e. the absorption lines of the approaching star appear 
stronger (Conti et al.[T9"75 Seggewiss & de Groot l 19761 Lloyd 
Evans I1979> . although the reverse effect is seen in the UV 
(Stickland et al. 1993 ). Though the origin of this effect is as yet 
not established, a commonly proposed scenario involves the ex- 
istence of an interaction process within the binary system (Gies 
et al. lT997l . 

Analysing UV resonance line profiles of HD 159176 as ob- 
served with IUE, Pachoulakis ( 1996) derived a mass-loss rate 
of about 3 10~ 6 M Q yr _1 for each star (note that this value is 
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a factor five larger than the one derived by Howarth & Prinja 
I1989> . The stellar winds of the components of HD 159176 are 
therefore probably sufficiently energetic to interact. The inter- 
action region is expected to be located roughly mid-way be- 
tween the stars and it prevents the wind of each star from de- 
ploying into the direction towards the other star. This situation 
has an impact on the UV resonance line profiles and the UV 
light curve (Pachoulakis lT996l Pfeiffer et al. [T997l . Pfeiffer et 
al. analysed the variations of the UV resonance lines and con- 
cluded that there exists a source of extra emission between the 
stars. These authors associated this extra emission with reso- 
nant scattering of photospheric light by the material inside a 
colliding wind region. According to their analysis, this shock 
region is slightly wrapped around the secondary star. 

HD 159176 was detected as a rather bright X-ray 
source with EINSTEIN (0.157 ± 0.008 ctss" 1 with the IPC, 
Chlebowski et al. [T9891 and ROSAT (0.291 ± 0.034CISS- 1 
during the All Sky Survey, Berghoferet al. 1996 1. Chlebowski 
& Garmany ( 1991 ) suggested that the excess X-ray emission 
observed in many O-type binaries compared to the expected 
intrinsic contribution of the individual components is produced 
by the collision of the stellar winds (see also e.g. Stevens et 
al. [1992 1. Therefore, it seems likely that at least part of the 
X-ray flux of HD 159176 may originate in the wind interaction 
region. In contrast with this picture, Pfeiffer et al. ( 1997 > 
suggest that the bulk of the X-ray emission arises primarily 
from the intrinsic emission of the individual components, 
rather than from a colliding wind interaction. 

To clarify this situation, we obtained an AOl XMM-Newton 
observation of HD 159176. Since the system is X-ray bright, 
it is well suited to investigate the X-ray properties of a short- 
period early-type binary. 

Table 1. Relevant parameters of the HD 159176 binary sys- 
tem adopted throughout this paper unless otherwise stated. The 
numbers are taken from Pachoulakis ( 1996 P96) and Diplas & 
Savage (TT994I DS). 



Parameter 


Prim. 


Sec. 


Ref. 


a sin i (Rq ) 




28.9 


P96 


M(M ) 


31.9 


31.6 


P96 


in 




~ 50 


P96 


Teff (K) 


42500 


35000 


P96 


R* (Re) 


9.8 


9.3 


P96 


Voo (kms -1 ) 




2850 


P96 


M(M yr" 1 ) 


3.2 10" 


6 2.6 HT 6 


P96 


log N h ,ism (cm -2 ) 




21.23 


DS 



2. Observations 

The observation of HD 159176 with XMM-Newton (Jansen et 
al. 2001 ) took place during revolution 229, in March 2001 
(JD 245 1977.903- 2451978.338). According to the ephemeris 
given by Stickland et al. dl993> . the whole observation covered 



phases between 0.53 and 0.66 (phase 0.0 corresponding to the 
maximum radial velocity of the primary component). 

The total observation was split into three exposures. During 
the first exposure (Obs. ID 0001730401), the satellite was still 
inside the radiation belts and only the RGS instruments were 
used for on-target observations. The total on-target EPIC ex- 
posure time of the other two exposures (Obs. IDs 0001730201 
and 0001730301) is about 26 ks. 

The three EPIC instruments were operated in the full frame 
mode (Turner et al. l200Tl Striider et al. 2001 1. All three EPIC 
instruments used the thick filter to reject optical light. The two 
RGS instruments were operated in the default spectroscopy 
mode (den Herder et al. 2 001t . 

We used version 5.2 of the XMM-Newton Science Analysis 
System (SAS) to reduce the raw EPIC data. For each EPIC cam- 
era, the observation data files (ODFs) of the two exposures 
were merged by the XMM-Newton Science Operation Center 
(SOC) into a single set of observation data files. In order to 
process the EPIC-MOS files with the 'emproc' pipeline chain, 
we had to manually correct some keywords in the ODF sum- 
mary file generated by the 'odfingest' task. For the pn camera, 
the merged ODF had to be reprocessed by the XMM-Newton 
SOC to generate an event list. 

The only events considered in our analysis were those with 
pattern - 12 for EPIC-MOS and for EPIC-pn respectively 
(Turner et al. 2001 ). We found no event with pattern 26-31, 
so that we conclude that no significant pile-up affects our data. 

In this paper, we discuss only the X-ray data for 
HD 159176. The other X-ray sources found in the EPIC field 
of view and associated with the open cluster NGC6383 are 
discussed in a separate paper (Rauw et al. 2003). The response 
matrices used for the EPIC spectra of HD 159176 were those 
provided by the SOC. No difference was found between spec- 
tra obtained with these response matrices and those obtained 
with the SAS-generated ones. 

For the RGS data, separate ODFs were provided for each 
of the three observations. The three data sets yield exposure 
times of respectively about 10, 16 and 9ks. The RGS data were 
processed with version 5.3 of the SAS. The raw data were run 
through the 'rgsproc' pipeline meta-task. Appropriate response 
matrices were generated using the 'rgsrmfgen' SAS task. 

3. Dealing with high background level episodes 

We extracted a light curve at very high energies (Pulse 
Invariant (PI) channel numbers > 10000) that revealed a high 
level of soft proton background. The soft protons responsible 
for these flares are thought to be accelerated by magneto- 
spheric reconnection events unrelated to solar flares (Lumb 
2002). Unfortunately this high-level background affects a 
significant fraction of the exposures. Indeed, nearly one half of 
the total EPIC observation is affected. For this reason, one has 
to check whether a bad time interval rejection is necessary, or 
whether the background correction of spectra could take these 
undesired events properly into account. As a convincing test, 
we generated EPIC spectra on the basis of the whole event 
lists (including events occurring during the flare). On the other 
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Table 2. Performed and useful exposure times for all XMM- 
Newton instruments. The useful RGS exposure time is con- 
siderably reduced due to a strong flare of soft protons. The 
fourth column yields the background-corrected count rate for 
the source, between 0.4 and 10 keV for EPIC, and between 
0.3 and 2.5 keV for RGS. In the case of the RGS instruments, 
only exposure 0001730201 is considered. The error bars on the 
count rate represent the ± 1 a standard deviation. 



Instrument 


Performed 


Effective 


Background-corrected 




duration 


exposure 


count rate 




(s) 


(s) 


(counts - ) 


MOS1 


25 901 


25 435 


0.512 ±0.005 


MOS2 


25 906 


25 473 


0.532 ± 0.005 


pn 


22759 


20402 


1.303 ±0.009 


RGS1 


16495 


7490 


0.085 ± 0.006 


RGS 2 


16535 


6838 


0.101 ±0.006 




1.0054x10° 

Time (s) 



Fig. 1. Light curve for events with PI > 10000 obtained for 
the EPIC-MOS2 field and covering the two exposures where 
the EPIC instruments were operating. A time bin of 100 s was 
adopted. 



hand, we filtered the event lists for time intervals with high- 
energy (PI > 10000) count rates exceeding thresholds of 0.20 
and l.lOcts s _1 , for the MOS and pn cameras respectively, in 
order to reject the flares and obtain filtered spectra. Of course, 
the spectra obtained with the complete event list have a better 
signal-to-noise ratio than the spectra obtained from the filtered 
one because of the longer exposure time (about 25 ks for the 
whole exposure, instead of about 13 ks for the filtered data). 
Following detailed comparison of the results from the two 
approaches, we did not find any significant difference between 
the fitting results (see Section^ relevant to both cases. For this 
reason, we decided not to reject the flare intervals from our 
EPIC data set, and the effective exposure times (see Table|2j 
are thus close to the performed observation duration. 

Timing analysis of the flare events was also performed to 
evaluate their impact on the lower energy (PI<10000) light 
curves. We extracted EPIC-pn events with boxes covering the 
12 CCDs, but excluding the sources appearing in the field. 
Light curves were extracted, between PI 400 and 10000, with 
a PI width of 400 units and time bins ranging from 50 to 
1000 s. These light curves show that the mean level of flare 
contribution to the background is only about twice the mean 
level of the proton unaffected background. Fig.|2] shows the 
count rate versus PI, where each point stands for a 400 PI 
channel interval, for the time interval affected by the main 
flare. The count rate of this background decreases by about 
a factor 10 from PI 400 to 10000. The bump near PI 8000 is 
due to fluorescent emission lines produced by interaction of 
charged particles with the camera body material (Lumb 2002 1. 
Comparing the overall background level (including the flare) 
to the count rate in the source region over PI G [400:800], i.e. 
over the energy range with the highest background level, we 
find that, at its peak level, the total background contributes 
only for about 0.9 % of the corrected source count rate. The 
same calculations performed in other energy bands reveal 
that the background contribution below 5 keV is always at 
least 4 times lower than the source count rate. This result 
confirms that the rejection of the soft proton time intervals is 
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Fig. 2. EPIC-pn count rates versus PI for the global background 
(including the soft proton flares) for a ~ 350 arcmin 2 area 
(about a half of the whole detector plane area). The background 
is much stronger in low energy bands than at higher energies. 
The bump seen near PI = 8000 is due to fluorescent emission 
lines produced by charged particles interacting with the cam- 
era body material. The error bars represent the ± 1 a standard 
deviation. 



irrelevant for the analysis of HD 159176, and that the flare is 
properly taken into account by the background subtraction step. 

However, this conclusion does not hold for RGS data. 
Indeed, the comparison of spectra obtained with and without 
the flare epochs revealed significant differences. As a conse- 
quence, low-background good time intervals were selected 
by application of count thresholds to the high energy light 
curves. These thresholds were evaluated, on the basis of a 
visual inspection of the light curves. The third (9ks) exposure 
was almost completely contaminated by high background, 
whereas the first (10 ks) observation gave only about 3ks of 
useful time. Only the 16ks exposure could be used, although 
more than half of it had also to be rejected. For this reason, the 
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effective RGS exposure time is severely reduced, as shown in 
TableEl 



4. The EPIC spectrum of HD 159176 

The source spectrum was extracted within a circular region 
of 60arcsec radius. The intersection with a circular region 
of radius 15arcsec centered on a faint point source located 
at RA=17:34:40 and DEC=-32:34:33 (Equinox J2000) was 
rejected. The background was selected within an annular 
region centered on HD 159176, with an outer radius 1.414 
times larger than for the source region. Again, the intersections 
with 15 arcsec circles centered on two point sources (at 
RA=17:34:48, DEC=-32:34:24 and RA=17:34:48, DEC=- 
32:35:21) were rejected. Consequently, source and background 
regions extend nearly over the same area. We also tried to 
select alternative background regions with boxes distributed 
near the source region, and avoiding the numerous point 
sources of the field. This selection did not give a significant 
difference in the background corrected spectrum. In the case 
of the pn camera, an additional difficulty arises from the gaps 
between the CCDs of the pn detector. We excluded the gaps 
manually from the source and background regions. 
The additional figure (jpeg format) shows the EPIC-MOS1 
field with the selected regions for the source and the back- 
ground in the annular case. The small circular regions are 
used for the exclusion of the point sources located inside these 
source and background regions. 

The emission lines of the RGS spectrum (see the next 
section) indicate that at least the low energy X-ray emission of 
HD 159176 is dominated by a thermal plasma emission. On 
the other hand, the shock region in a colliding wind binary 
might accelerate free electrons up to relativistic velocities and 
in combination with the intense UV radiation fields of the 
two stars, these particles could produce a hard non-thermal 
X-ray emission through inverse Compton scattering (see Chen 
& White lT99"Tl Eichler & UsovT993). Therefore, we tested 
various models combining optically thin thermal plasma 
models with a non-thermal power law component. In the 
framework of a model where the X-ray emission arises from 
shocks (either in the wind interaction or in instabilities of 
the individual winds), one should observe a range of plasma 
temperatures corresponding to the distribution of the pre-shock 
wind velocities. Therefore, one expects a priori that more 
than one thermal component is required to fit the observed 
spectrum. So, the spectra have been fitted either with two 
temperature mekal thermal plasma models (Mewe et al. 1985 
Kaastra 1992 1 or with models containing a single temperature 
mekal component along with a non-thermal power law 
component. For each thermal component, solar metallicity was 
assumed. 

The spectra were rebinned to achieve a minimum of 9 and 25 
counts per energy bin for the MOS and pn data respectively 
and then fitted with the XSPEC software (version 11.1.0). For 
each model, we used an interstellar neutral hydrogen column 
density of N H = 0.17 10 22 cm" 2 (Diplas & Savage IT9941 . and 




Energy (keV) 

Fig. 3. EPIC-MOS1 spectrum of HD 159176 fitted with the 
two-temperature mekal model between 0.4 and 5.0 keV. 



a second absorption component to account for circumstellar 
absorption. The results are listed in Table|3] We note that the 
best-fitting power law components turn out to be rather steep 
with photon indices (V) above 3.8. This situation is very much 
reminiscent of the results obtained for the 04 V star 9 Sgr 
(see Rauw et al. I2002b> . Figure [3] illustrates the spectrum 
of HD 159176 as observed with the MOS1 instrument. Note 
that including different absorption columns for the individual 
model components does not improve the quality of the fits. In 
fact, XSPEC fails to fit at least one of the column densities. 

Another model was fitted with some success to the 
EPIC-pn data. This model consists of three components: 2 
thermal mekal components and a non-thermal power law 
component (see Sect.|8] for a possible interpretation of this 
component). This fit yields a better x^, albeit the value 
obtained for the photon index is once more rather large. The 
results are summarized in Table|4] for the two background 
regions selected. One can see that both background techniques 
yield very similar results. This comforts us with the idea that 
background selection is not a critical issue for this source, and 
that the annular region with parasite sources excluded provides 
a satisfactory description of the actual background. Figure|4] 
displays the pn spectrum along with the best-fit corresponding 
to this latter model. We note that this three component model 
yields a significantly better fit to the EPIC-pn spectrum above 
2 keV than the two component models. 

The main motivation behind these fits is to provide an 
overall description of the observed spectral distribution that 
can then be used to derive fluxes and to check for consis- 
tency with previous observations. With the parameter values 
obtained for these fits, we can determine the flux emitted by 
HD 159176 between 0.4 and 10.0 keV. For this purpose we 
adopt the best fitting model obtained for EPIC-pn, i.e. the three 
component model whose parameters are shown in Table|4] The 
flux corrected for interstellar absorption only is 5.52 10 -12 
ergcm~ 2 s -1 . If we adopt a distance to NGC6383 of 1.5 kpc, 
we obtain a luminosity of 1.48 10 33 ergs -1 . This value can 
be compared to X-ray luminosities obtained with EINSTEIN 
(Chlebowski et al. lT9S9l and ROSAT (Berghofer et al. H996l . 
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Table 3. Fitted model parameters for EPIC spectra of HD 159176. The upper part of the table concerns 
the wabsisM*wabswiND* (mekal + mekal) model, whilst the bottom part summarizes the results for the 
wabsisM*wabswiND* (mekal + power) model. The Nh given in the table refers to the circumstellar (i.e. wind) column 
assumed to be neutral. The interstellar column density is frozen at 0.17 1 22 cm -2 for all models. The Norm parameter of the 
mekal component is defined as (10~ 14 /(4 it D 2 )) J n e nu dV, where D, n c and are respectively the distance to the source 
(in cm), and the electron and hydrogen densities (in cm -3 ), whereas for the power law this parameter corresponds to the photon 
flux at 1 keV. The error bars represent the 1 a confidence interval. 







kT 1 


Norm 1 


kT 2 


Norm^ 


v 2 




(10 22 cm" 2 ) 


(keV) 




(keV) 




d.o.f. 


MOS1 


0.41 


0.22 


2.04 10~ 2 


1.05 


1.51 10~ 3 


2.19 




±0.03 


±0.04 


±0.38 10" 2 


± 0.02 


± 0.09 10~ 3 


197 


MOS2 


0.39 


0.21 


2.11 10" 2 


1.00 


1.58 10 -3 


2.65 




±0.03 


±0.04 


± 0.40 10~ 2 


± 0.02 


± 0.09 10~ 3 


198 


pn 


0.35 


0.18 


1.88 10~ 2 


0.59 


2.98 10 -3 


1.87 




±0.02 


±0.03 


± 0.26 10~ 2 


± 0.01 


±0.17 10 -3 


321 




N H 


kT 


Norm 1 


r 


Norm 2 


xi 




(10 22 cm -2 ) 


(keV) 








d.o.f. 


MOS1 


0.26 


0.30 


4.69 10 -3 


4.01 


1.87 10 -3 


1.65 




±0.02 


±0.01 


±0.84 10~ 3 


±0.09 


±0.18 10~ 3 


197 


MOS2 


0.24 


0.29 


5.08 10~ 3 


3.80 


1.69 10" 3 


1.99 




±0.02 


±0.01 


±0.89 10" 3 


± 0.08 


±0.16 10~ 3 


198 


pn 


0.28 


0.28 


6.72 10 -3 


4.09 


1.67 10 -3 


1.70 




±0.01 


±0.01 


±0.80 10~ 3 


±0.08 


±0.14 10 -3 


321 



Table 4. Fitted model parameters for the EPIC-pn spectrum of HD 159176 for two different background estimation techniques. 
The adopted model is wabsisM*wabswiND* (mekal + mekal + power ). The column densities given in the table refer 
to the circumstellar material. The interstellar column density is frozen at 0.17 1 22 cm -2 . The error bars represent the 1 a 
confidence interval. 





N H 


kT 1 


Norm 1 


kT 2 


Norm 2 


r 


Norm 3 


xi 




(10 22 cm' 2 ) 


(keV) 




(keV) 








d.o.f. 


Annulus 


0.20 


0.20 


4.39 10~ 3 


0.58 


1.79 10" 3 


3.50 


7.35 10"" 


1.11 




±0.02 


±0.01 


± 1.05 10" 3 


±0.01 


±0.19 10" 3 


± 0.12 


± 1.25 10~ 4 


319 


Boxes 


0.18 


0.21 


3.66 10~ 3 


0.59 


1.72 10" 3 


3.43 


8.12 10~ 4 


1.09 




±0.01 


±0.01 


± 0.83 10" 3 


±0.01 


±0.18 10" 3 


±0.10 


± 1.13 10~ 4 


319 



In the case of EINSTEIN, Chlebowski et al. converted the IPC 
count rates between 0.2 and 3.5 keV into fluxes, assuming 
a 0.5 keV thermal bremsstrahlung model and an absorption 
column (0.209 10 22 cm~ 2 ) calculated on the basis of the color 
excess E(B — V) determined through a relation given by 
Mihalas & Binney ( 198 1 ) for an intrinsic color index (B — V)q 
of —0.309 and an observed (B — V) color index of ±0.04. 
For ROSAT, Berghofer et al. converted the PSPC count rate 
(between 0.1 and 2.4 keV) into a flux assuming a 0.73 keV 
Raymond & Smith dl977> optically thin plasma model. The 
same absorption column as for EINSTEIN was used. The X-ray 
luminosities scaled to the 1.5 kpc distance are shown in the 
second column of Table|5]for EINSTEIN, ROSAT, and XMM. 
The luminosities derived from the different satellites are quite 
different. These discrepancies cannot be explained by the 
differences in the sensitivity ranges of the various instruments. 
Instead, it seems more likely that the model assumptions 
used to convert the EINSTEIN-IPC and ROSAT-PSPC count 
rates into luminosities are responsible for the discrepant 




Energy (keV) 

Fig. 4. EPIC-pn spectrum of HD 159176 fitted with the two 
mekal components + power law model between 0.4 and 
5.0 keV. The three components are also individually displayed. 
The power law begins to dominate the spectrum at energies 
above 2 keV. 
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luminosities. To check this, we have used XSPEC to simulate 
EINSTEIN-IPC and ROSAT-PSPC spectra of HD 159176. The 
three component model derived from the fit of the EPIC-pn 
data (see Table 4) was folded through the response matrices of 
the IPC and PSPC instruments. The corresponding theoretical 
count rates, 0.138 cts s _1 for the IPC and 0.287 cts s _1 for the 
PSPC, are in very good agreement with the observed values 
(see Table 5). Therefore, we conclude that the archive data do 
not reveal evidence for long-term variability of the X-ray flux 
of HD 159176. 

Since the orbit of HD 159176 does not have a significant 
eccentricity, it is expected that a possible orbital modulation 
of the X-ray emission from a colliding wind interaction zone 
would be due entirely to the changing circumstellar column 
density along the line of sight. However, given the rather low 
inclination (i ~ 50°, Table 1), we expect only a marginal 
orbital modulation of the X-ray flux at energies above 0.4 keV. 
At this stage, it is worth recalling that our XMM observation 
of HD 159176 was obtained during orbital phases 0.53-0.66, 
i.e. shortly after quadrature. Therefore, we expect no strong 
occultation effect of the X-ray emission from either component 
or from the wind interaction region: our value of the luminosity 
should thus reflect the typical X-ray emission of HD 159176. 

We use the total bolometric luminosity of HD 159176 as 
derived from the parameters in Tabled i.e. 1.49 10 39 ergs -1 , 
to compute the Lx/Lboi ratios given in Table|5] On the other 
hand, we estimate the X-ray luminosity of the individual 
components by means of the empirical relation of Berghofer 



et al. J1997i . The ratio of the sum of the expected X-ray 
luminosities and the total bolometric luminosity amounts to 
1.48 10~ 7 , which is about a factor 7 inferior to the XMM value 
listed in Table[5] suggesting that there exists indeed a moderate 
excess X-ray emission in HD 159176. 



Table 5. X-ray luminosities and luminosity ratios from 
EINSTEIN, ROSAT and XMM observations of HD 159176. 
Luminosities are calculated for a 1 .5 kpc distance. The Obs. CR 
column gives the count rates from Chlebowski et al. ( 1989 1 and 
Berghofer et al. (H9961 respectively for EINSTEIN and ROSAT, 
as well as our value for EPIC-pn. The Th. CR column provides 
the theoretical count rates calculated on the basis of the individ- 
ual instrument response matrices for the 3 component model of 
Table[4](see text). 



Satellite 


Lx 


Lx/Lboi 


Obs. CR 


Th. CR 




(ergs -1 ) 




(ctss -1 ) 


(ctss -1 ) 


EINSTEIN 


4.33 10 33 


2.89 10~ B 


0.157 


0.138 


ROSAT 


8.43 10 32 


5.64 10 -7 


0.291 


0.287 


XMM 


1.48 10 33 


9.93 1(T 7 


1.303 


1.328 



5. The RGS spectrum of HD 1 591 76 

Because of the strong soft proton activity and since it turned 
out to be impossible to merge event lists from separate RGS 




1 5 20 
Wavelength (A) 



Fig. 5. Combined fit of RGS 1 and RGS2 spectra between 8 and 
27 A for a two-temperature mekal model. The strongest lines 
are labelled. 



exposures, we had to restrict our analysis to the results of the 
pipeline processing of the 16ks exposure. We extracted first 
order spectra of HD 159176 via the 'rgsspectrum' SAS task. 
Figure[5] shows RGS1 and RGS2 spectra. The most prominent 
features in the RGS spectrum are the Lya lines of Ne X and 
O VIII, as well as the He-like triplet of O VII and some strong 
Fe XVII lines. The Ne IX He-like triplet is also present, but 
the poor quality of our data as well as the blend with many 
iron lines from various ionization stages hamper the analysis 
of this part of the spectrum. The temperatures of maximum 
emissivity of the lines seen in the RGS spectra of HD 159176 
span a range from - 2 to - 6 10 6 K (i.e. kT ~0.17 - 0.52 keV, 
see the APED database, Smith & Brickhouse 2000). 

Models were fitted to the RGS spectra. The model includ- 
ing the two-temperature mekal component yields the best 
fits. The combined fit of RGS1 and RGS2 spectra by a two- 
temperature mekal thermal model is also illustrated in Fig. [5] 
The best fitting parameter values are given in Table|6] With 
reasonable xt values, we obtain temperatures (about 0.16 and 
0.60 keV) consistent with the fit parameters of the EPIC pn 
spectrum. It is interesting to note that these temperatures cor- 
respond to the extreme values of the range inferred hereabove 
from the maximum emissivities of the lines identified in the 
spectrum. The interstellar column density appears also com- 
parable within the error bars to the EPIC results for the same 
two-temperature model. 

The agreement between the RGS and EPIC pn parameters 
suggests that we can attempt a simultaneous fit of the spec- 
tra from different instruments. For this purpose, we used the 
two-temperature mekal model. The results are summarized in 
Table|6] First, we fitted the three EPIC spectra together. Next, 
we also included the RGS spectra. 

We see that, in the first case, best fitting parameter values 
are very close to the results obtained hereabove for the MOS 
spectra (which are somewhat different from the best fitting 
parameters of the pn spectrum). On the other hand, including 
also the RGS spectra in the fit, we obtain a value of the 
second temperature which is also in good agreement with 
that obtained for the MOS spectra fitted individually. Finally, 
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Table 6. Model parameters for a two-temperature mekal model fitted to the RGS1 and RGS2 spectra of HD 159176, to 
the three EPIC spectra considered simultaneously and finally to all the XMM data taken together. The adopted model is 
wabsisM*wabswiND* (mekal + mekal) . The column density given in the Table refers to the local (i.e. circumstellar) 
absorption. The interstellar column density is frozen at 0. 17 1 22 cm -2 . The error bars represent the 1 a confidence interval. 





N H 


kT 1 


Norm 1 


kT 2 


Norm 2 


xi 




(10 22 cm -2 ) 


(keV) 




(keV) 




d.o.f. 


RGS1 


0.32 


0.16 


1.66 10~ 2 


0.62 


3.96 10-' 1 


1.36 




±0.11 


±0.01 


± 1.70 10~ 2 


±0.04 


± 1.81 10 -3 


61 


RGS2 


0.27 


0.16 


1.37 10~ 2 


0.59 


2.67 10 -3 


1.45 




±0.09 


±0.01 


± 0.26 10~ 2 


±0.03 


± 1.05 10" 3 


71 


EPIC 


0.40 


0.21 


2.19 10 -2 


0.97 


1.48 10 -3 


2.41 




±0.14 


±0.02 


± 0.23 10~ 2 


±0.01 


± 0.05 10~ 3 


726 


EPIC 


0.40 


0.21 


2.22 10~ 2 


0.96 


1.54 10" 3 


2.28 


+ RGS 


±0.01 


±0.02 


± 0.22 10~ 2 


±0.01 


± 0.05 10~ 3 


862 



we note that fitting the EPIC + RGS data with a differential 
emission measure model (i.e. the c6pmekl model in XSPEC, 
see Singh et al. 1996 1 yields a differential emission measure 
with a very broad maximum spanning from about 0.15 to 
1.2keV, again in reasonable agreement with the temperature 
ranges inferred from the RGS line spectrum. 
The best fitting temperatures of the X-ray spectra of 
HD 159176 are rather low. In fact, these temperatures are 
roughly typical for presumably single O-stars. On the other 
hand, for an O + O colliding wind binary, one expects a priori 
to find harder X-ray emission from the wind interaction region. 
For instance, in the spectrum of HD 93403, Rauw et al. ( 2002a ) 
found a hot thermal component with kT ~ 0.92 - 1.25 keV, 
which is significantly hotter than the temperature of the harder 
component in the EPIC and RGS spectra of HD 159176. 
However, the smaller orbital separation of HD 159176 (as 
compared to HD 93403) prevents the stellar winds of its 
components from reaching their terminal velocities before they 
collide, leading to lower temperatures in the wind interaction 
region. 

Given the rather poor signal-to-noise ratio of our RGS data, 
nothing reliable can be said about the morphology of the emis- 
sion lines. For instance, the strongest line of the spectrum, the 
O VIII Ly a line at 18.97 A, is illustrated in Fig.|6] This line ap- 
pears to be significantly broadened. In fact, the line's FWHM 
estimated from Fig.[6]is about 2500km s _1 , while the FWHMs 
of the instrumental profiles at this wavelength are about 0.06 
and 0.07A (i.e. 950 and 1100 kms" 1 ) for RGS1 and RGS2 
respectively (den Herder et al. 2001 1. The line appears there- 
fore significantly larger than what would be expected for a 
coronal emission model (see Owocki & Cohen 2001 ). We also 
note that the line profile is neither gaussian nor lorentzian. 
Theoretical models of line emission from either single O-type 
stars or colliding wind binaries (e.g. Owocki & Cohen 12001 1 
Feldmeier et al. 2003 Henley et al. |2003 1 predict rather flat- 
topped profiles with asymmetries depending on the amount of 
wind absorption and, in the case of binary systems, on the or- 
bital phase. The detection of an orbital modulation of the X- 
ray line profiles such as predicted by Henley et al. would pro- 



vide a clear evidence that the lines are formed in the wind 
collision zone. However, since we have only a single XMM- 
Newton observation of HD 159176 at our disposal, we cannot 
make any statement about the existence or not of such a mod- 
ulation. Unfortunately, our current data set does not have the 
quality needed to perform a detailed comparison between the 
observed lines and the predictions from the various models. All 
we can state is that the O VIII Ly a line is slightly blue-shifted 
by about 300 - 600 km s , though the highest peak in the pro- 
file is found near the rest wavelength. Within the error bars, 
these features are in reasonable agreement with some of the 
line profiles predicted by Owocki & Cohen ( 2001 ) for a single 
O star wind. On the other hand, the agreement with line pro- 
files formed in a wind interaction zone for a wind momentum 
ratio of r\ ~ 1 (see below) and an orbital phase of 0.53 - 0.66 
(corresponding to 8 6 [98°, 130°] in the formalism of Henley 
et al. 2003) is slightly less good. At these phases, the colliding 
wind models predict the highest peak to occur with a signifi- 
cant blue-shift (Henley et al.). We caution however that the line 
profiles presented by Henley et al. were computed assuming an 
adiabatic wind interaction, while the wind interaction zone in 
HD 159176 is most likely radiative (see below). 

Another noteworthy feature of the RGS spectra of 
HD 159176 is the fact that the forbidden (/) component 
of the Ovu He-like triplet is very weak. Inspection of the 
spectrum reveals no significant excess of counts above the 
level of the pseudo-continuum at the wavelength of the / 
line, whereas the intercombination (i) and resonance (r) 
lines are clearly detected at a level of about 7 x 10~ 5 and 
3 x 10~ 5 photons cm~ 2 s _1 respectively. Unfortunately, the 
limited quality of our data does not allow us to perform 
sophisticated fits of the triplet. Nevertheless, we have fitted the 
triplet by three gaussians plus a constant (to mimic the pseudo- 
continuum) where we requested that the three gaussians should 
have the same Doppler shift and the same width. This fit yields 
an / ji ratio of 0.08 which can probably be considered as a safe 
upper limit. A low / /i ratio is not unexpected. In fact, in the 
winds of hot early-type stars, the intense UV radiation induces 
a strong coupling between the upper level of the forbidden line 
and the upper level of the intercombination line, leading to 
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Fig. 6. The Ovill Ly a line at 18.97 A as seen with the RGS1 
(dotted line) and RGS2 (solid line) instrument. 



an enhancement of the i and a decrease of the / component 
(Porquet et al. 1200 H . For the O VII He-like triplet, Porquet 
et al. calculate an f/i ratio of less than 3 x 10~ 3 for a star 
with a radiation temperature > 30 000 K and for a dilution 
factor w > 0.1. In HD 159176, this mechanism can act on the 
intrinsic X-ray emission of the individual binary components 
as well as on the X-ray emission from the colliding wind 
region. Indeed, assuming that the wind-interaction zone lies 
about half-way between the two stars, the UV radiation field of 
each star is diluted by w = 0.15 at the location of the colliding 
wind region. The enhanced plasma density in the colliding 
wind region could further reduce the / ji ratio (Porquet et al. 
2001 1, though this is likely to be only a second order effect in 
the case of a hot, close binary such as HD 159176. 

Finally, a further comparison of the RGS spectrum (Fig.[5J 
with RGS spectra obtained for stars such as 9 Sgr (Rauw 
et al. l2002b> or (Pup (Kahn et al. l2001i shows that the 
HD 159176 spectrum does not display the N VII A 24.8 A line. 
To find out whether this is due to a reduced N abundance or not, 
we fitted simultaneously RGS 1 and RGS2 spectra with variable 
abundance mekal models, allowing the nitrogen abundance to 
vary (with other abundances frozen at the solar value). The re- 
sulting fit suggests a slightly (2.85 ± 1.26) enhanced N abun- 
dance (with respect to solar abundance) while we recover the 
temperatures for the RGS fits within the error bars. Moreover, 
the N VII A 24.8 A line appears clearly in the theoretical RGS 
spectrum of HD 159176 when the absorption columns are set 
to zero, albeit this line nearly completely disappears when the 
absorption components are set to the values derived in our fit. 
Therefore, the absence of the N VII line in our observed spec- 
trum is due to the strong absorption that photons at this rather 
low energy suffer when they encounter the circumstellar or in- 
terstellar medium. Also, the relative intensities of oxygen lines 
compared to other elements seem quite higher in the spectrum 
of HD 159176. At this stage, it is worth recalling that CPup 1S 
an evolved object and the line intensities in its RGS spectrum 
reflect the signature of the CNO process (Kahn et al. 2001 1. On 



the other hand, we do not expect to see such an effect in the 
spectrum of a relatively unevolved system such as HD 159176. 



6. Short-term variability of the X-ray flux 

Hydrodynamic models of colliding wind binaries predict the 
existence of various kinds of dynamical instabilities that should 
affect the wind interaction region (Stevens et al. 1992, Pittard 
& Stevens fl 9971 1 and manifest themselves through a rapid vari- 
ability of the emerging X-ray flux in addition to the much 
slower orbital phase-locked modulation. Theoretical consider- 
ations (Stevens et al. 1992 1 indicate that the nature of these in- 
stabilities depends on the efficiency of radiative cooling in the 
shocked region which is a function of the orbital separation and 
of the wind properties of the components of the binary system. 
Previous X-ray satellites lacked the required sensitivity to look 
for the signature of these variations. Owing to the large collect- 
ing area of the X-ray mirrors onboard XMM-Newton and since 
HD 159176 is rather X-ray bright, we can use our present data 
set to seek for the signature of these instabilities. 



1 .0053* 1 O" 



1.0054X 1 3 



1 .0055x 10 3 



0.6 - i 

- i * 



i i i I * i I i J 1 



i ■ 



1 .0053* 1 O a 



1.0054X 1 O a 



1.0055X 10 3 



* * i 



• + 



1 .0053* 1 O a 



1 .0054* 1 O e 



1 .0055x 10 3 



Fig. 7. Count rates evaluated in the 0.4-10.0 keV energy band 
versus time for MOS 1 (upper panel), MOS2 (middle panel) and 
pn (bottom panel). The time bin used is 1000 s for all instru- 
ments. The error bars represent the ± 1 a standard deviation. 
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To investigate the putative short-term variability of 
HD 159176, we divided the effective exposure into short time 
intervals. In the case of the MOS instruments, we used 2.1 ks 
time intervals. We obtained 9 and 10 equal time intervals 
respectively for MOS1 and MOS2. Because the good time 
intervals (GTIs) of the different instruments are not strictly 
identical, the time bins of the time series of the two instruments 
do not match exactly. A spectrum was extracted for each time 
bin through the same procedure as described in Section^ Each 
spectrum was rebinned, to get at least 4 counts per energy 
channel, and fitted with a two temperature mekal model. 

Some parameters seem indeed to change with time but 
no correlation is found between MOS1 and MOS2. In each 
case, the lower temperature does not vary significantly while 
the higher temperature component displays larger variations. 
No systematic trend is found for these variations. The same 
study was carried out for EPIC-pn data using the two models 
described in Sect. 3 for time intervals of 2 ks. In this case, 
we do not observe any significant variation in the model 
parameters either. 

A more model-independent way to check for variability 
consists in studying the run of count rates for particular time 
bins, over different energy ranges: see for instance Fig. Q for 
count rates between 0.4 and 10.0 keV within 1 ks time bins. 
This figure does not reveal any significant correlation between 
the different EPIC instruments. The same conclusion applies 
also for narrower energy ranges selected within this same in- 
terval. Moreover, we see that count rates in all energy bands 
are nearly constant within 2 a for all instruments. 

Finally, we extracted light curves for the three EPIC instru- 
ments with smaller time bins of 200 s and 100 s, within spa- 
tial regions corresponding to the source and the background 
areas selected for spectra extraction (Sect.[4}. Because of stan- 
dard GTIs, some of the time intervals were shorter than 200 s 
(or 100 s) and we therefore scaled the count rates to the effec- 
tive exposure time within each bin. The background corrected 
count rates display some fluctuations but do not reveal any cor- 
relation between the various instruments, x 2 tests performed 
on these corrected light curves confirm this lack of short time 
scale variability. 

7. Confrontation with theoretical models 

To complete our study, we can compare the overall properties 
of the observed spectrum of HD 159176 with the predictions of 
colliding wind models. 

7. 1. Analytical and hydrodynamical estimates of the 
X-ray luminosity 

A relevant quantity which should be introduced here is the 
cooling parameter: 

^cool «3 G?7 

x = — = — 

tcsc M_7 



where t coo i, t csc , v%, d^, and M_7 are respectively the cooling 
time, the escape time from the shock region near the line of 
the centers, the pre-shock velocity in units of 1000km s , the 
distance from the stellar center to the shock in units of 10 7 km, 
and the mass loss rate in units of 10~ 7 M Q yr _1 (see Stevens 
et al. H992i . Depending on the binary parameters, the gas in 
a colliding wind region could be either adiabatic (x>l) or 
radiative (\ < 1). While present day hydrodynamical models 
probably provide a safe estimate of the X-ray emission in 
adiabatic winds, their use in a highly radiative situation is 
more problematic. In fact, for values of \ lower than about 
2 or 3, the numerical simulations fail to spatially resolve the 
cooling length in the wind interaction region. Moreover, the 
thin-shell instabilities that appear in the simulation of these 
shocks introduce a mixing of hot and cool material that further 
complicates the prediction of the X-ray emission. Therefore, in 
the radiative case, it seems more appropriate to use analytical 
results to estimate the X-ray luminosity. 

Turning back to HD 159176, we can estimate the value 
of the relevant cooling parameter using the binary parameters 
given in Tabled In a close binary such as HD 159176, the stel- 
lar winds do not reach Voo before they collide. The pre-shock 
velocity should be in the range 850 - 2000 kms -1 . The upper 
boundary of this range is obtained assuming a 'standard' wind 
velocity law with (3 = 0.6 and v^ = 2850kms~ 1 ; the lower 
boundary corresponds to the value expected if we account for 
radiative inhibition effects for equal stars with a separation of 
about 40 R© (see Fig. 4 of Stevens & Pollock 119941 . These 
pre-shock velocity values are derived along the line connecting 
the centers of the stars, where the bulk of the X-ray emission 
is generated. 

Therefore, adopting the parameters of Tabled and a pre- 
shock velocity lower than 2000 km s~\ the resulting value of x 
is typical of radiative winds. The corresponding intrinsic X-ray 
luminosity can be estimated for each component of the binary 
through the following relation taken from Pittard & Stevens 
(2002 1 : 

£x,intr~0.5SMv 2 

where 3 accounts for a geometrical factor (1/6 for equal 
winds) and a small inefficiency factor (hence we used S 
= 0.1), and v is the wind speed at the contact surface. In 
the conditions described hereabove, Lx,intr is significantly 
larger than our observed value. For a pre-shock velocity of 
2000 kms -1 , Lx,intr (~ 8 1 35 ergs _1 ) is more than a factor 
400 too high. Note that the intrinsic luminosity needs to be 
corrected for stellar wind absorption before it can be compared 
to the observed luminosity. However, as we will see below, 
this correction hardly exceeds a factor of a few tens and is 
therefore not sufficient to account for such a huge discrepancy. 
If we reduce the pre-shock velocity to 850 km s _1 , the intrinsic 
luminosity becomes Lx,intr ~ 10 35 ergs -1 , which is still 
much too large. 

In order for our predicted luminosity to match the observed 
value of Lx, we need to reduce the mass-loss rates. As a 
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second step, we therefore adopt M = 610 _7 M Q yr~ 1 for 
both stars as suggested by Howarth & Prinja dl989> . With 
a pre-shock velocity of 2000 kms -1 , \ amounts to ~3.5, 
allowing us to derive the X-ray flux from a hydrodynamical 
simulation. For this purpose, we used the code of Pittard 
& Stevens ( 1997 > which is based on the VH-1 numerical 
code (see Blondin et al. H990i . This code uses a lagrangian 
piecewise parabolic method to solve the partial differential 
equations of hydrodynamics followed by a remap onto a 
fixed grid after each time step (for details see the references 
in Pittard & Stevens 1997 1. The simulations do not include 
radiative driving and assume instead constant wind speeds. 
A theoretical attenuated X-ray spectrum was generated from 
a 3D radiative transfer calculation on the hydrodynamical 
results, and we found Lx.att ~ 2 10 34 erg s _1 . Hence it seems 
that M must be reduced below the Howarth & Prinja d!989i 
value in order to account for the observed Lx. 

In order to explore a wider range of the wind parameter 
space, we decided to test an alternative approach that allows 
us to overcome some of the practical limitations encountered 
hereabove. This is described in the following section. 



7.2. Steady-state colliding wind model 

The main problem when modelling the X-ray emission from 
a close binary system like HD 159176 is that the wind-wind 
collision can be highly radiative. Since it is difficult, if not 
impossible, to resolve the width of such shocks in hydrody- 
namical calculations, the necessary information to compute the 
X-ray emission is lost. For this reason, Antokhin et al. (|2003 
in preparation) have developed a steady-state model able to 
account for the physics of radiative shocks in a way suitable to 
lead to a relevant discussion of colliding wind X-ray emission 
from close binary systems. 

In this approach, Antokhin et al. (2003 1 completely 
neglect instabilities and orbital rotation. Assuming a strong 
radiative cooling, both shocks generated in the primary and 
the secondary winds are thin. The pre-shock wind velocities 
are specified by a (3— law, and the shape of the contact 
surface is computed. The wind density and normal velocity 
are calculated at the exact position of the shock. The ID 
hydrodynamic equations are then solved for the case of a 
radiative post-shock flow. The X-ray spectrum is computed, 
and the process is repeated along the contact surface to obtain 
an integrated spectrum. The absorption of the X-ray radiation 
by the winds and by the shocks is determined using warm 
material opacity tables obtained with the CLOUDY v.94.00 
code (http://thunder.pa.uky.edu/cloudy/!. This model neglects 
any dissipative process and assumes that all the kinetic energy 
of the winds is radiated, thus giving an upper limit to the actual 
X-ray luminosity. For further details about this model, we refer 
the reader to Antokhin et al. (2003 1. 



2 r 




Fig. 8. Emerging X-ray luminosity between 0.5 and 10.0 keV 
as a function of the mass loss rate and the terminal velocity as 
computed for a distance of 1 .5 kpc. The parameter space covers 
values between 2.5 10~ 8 and 6.0 10 -6 M Q yr _1 for the mass 
loss rate, and between 1700 and 2450 kms -1 for the terminal 
velocity. 

7.2.1 . Equal winds 

Following this approach, a grid of models was generated and 
converted into a table model suitable to be used within XSPEC. 
As a first step, the assumption of equal winds was adopted. 
The parameters of the model are the mass loss rate and the 
terminal velocity. The parameter space covers M values from 
2.5 10~ 8 to 6.0 10~ 6 M Q yr _1 , and v^ values between 1700 
and 2450 kms -1 . The luminosities computed by the model 
were converted into fluxes assuming a distance of 1.5 kpc. 
The model normalization which scales as the inverse of the 
distance squared should be unity if all the model assumptions 
are correct and the distance is indeed equal to 1.5 kpc. Fig.|8] 
shows the emerging (i.e. after absorption by circumstellar 
material) luminosity between 0.5 and 10.0 keV versus the 
mass loss rate and the terminal velocity. This figure shows 
clearly that the general trend followed by the luminosity is 
to decrease as the velocity decreases. For M values lower 
than about 8.5 10~ 7 M Q yr _1 , Lx increases as the mass loss 
rate increases, but it decreases with increasing mass loss rate 
for higher M values. This turnover in the behaviour of the 
luminosity versus M is due to the increase of the optical depth 
of the warm wind material as the wind density increases. 

To take into account the fact that intrinsic emission of the 
individual winds is supposed to contribute to the observed 
spectrum, we used a mekal thermal plasma component 
in addition to the colliding wind component. Interstellar 
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absorption is accounted for by using an absorption column 
frozen to the value used in Sect.|4] The best fit of our EPIC-pn 
spectrum requires an additional power law component to fit 
the high energy tail already mentioned in the framework of the 
EPIC-pn spectral analysis in Sect.0] The physical motivation 
of this power law is discussed in Sect.[8j. Fig. [9] shows the fit 
obtained between 0.4 and 10.0 keV where the three individual 
components are displayed. The fluxes of the three components 
in the 0.4-10.0 keV band are 1.5 10 12 ergs" 1 , 3.8 10~ 12 
ergs -1 and 5.7 10~ 13 ergs -1 respectively for the mekal, 
the colliding wind and the power law components. This fit 
corresponds to logM = -6.67 ± 0.09, v^ = 1905 ± 50 
kms -1 , and a normalization parameter of 0.24 ± 0.04. In the 
best fit model, the mekal component has a temperature of 
0.29 ± 0.01 keV, and the power law has a photon index of 
2.50 ± 0.22. While this model provides an acceptable fit to the 
data - the reduced \l 1S 1-14 for 319 d.o.f. - we note that the 
best fit value of the normalization parameter is a factor 4 lower 
than expected. Because of the membership of HD 159176 to 
the open cluster NGC6383, it seems quite unlikely that the 
distance could be off by as much as a factor 2. Moreover, in 
order to explain the low normalization parameter, the actual 
distance of HD 159176 would have to be larger by a factor 2 
whereas an estimate of the distance on the basis of the expected 
absolute magnitude of an 07V + 07V binary rather suggests 
that the distance of HD 159176 might be less than 1.5 kpc 2 . 
If the normalization value is frozen to one, the model clearly 
overestimates the data mostly between 1 and 5 keV, leading 
to an unacceptable x 2 . Therefore, the low normalization value 
probably means that the model predicts, for the derived log M 
and Vqo, a luminosity which is too large by about a factor of 
4. A similar result is obtained for a simultaneous fit of MOS 1 
and MOS2 data with the same model. The only difference is a 
somewhat higher terminal velocity (Voo = 2190 ± 50 km s -1 ), 
compatible with the higher temperature obtained for mekal 
models described in Sect.|4] and a stronger disagreement of the 
normalization parameter which is a factor about 8 too low. 

The results of this subsection can be summarized as 
follows. Exploring a domain of the parameter space that 
includes the full range of reasonable values of M and Vqo for 
stars of spectral type 07 ± one subclass, we have found a 
model where the interplay between emission and circumstellar 
absorption allows to reproduce the observed spectral shape. 
However, this model as well as all other combinations of M 
and Vqo in this domain of the parameter space fail to reproduce 
the observed luminosity. The most likely explanation for this 
failure is that one or several of the model assumptions break 
down in the case of HD 159176. One possibility that we briefly 

1 We explored higher velocity values to try to fit this high energy 
tail with the colliding wind model. However, these models yield an 
increase of the intensity of the Fe K line at about 6.7 keV up to a level 
that is totally inconsistent with our data. 

2 With m v =5.64 (FitzGerald et al. U978> corrected for the redden- 
ing (Av=1.12) and Mv=^1.8 (Conti & Alschuler fT97H for a typical 
07V star, leading to Mv=-5.55 for an 07V + 07V system, we infer a 
distance of about 1.0 kpc only. With M v =-5.2 (Schmidt-Kaler 1982 1, 
the calculated distance is about 1.2 kpc. 
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Fig. 9. EPIC-pn spectrum of HD 159176 fitted with a mekal + 
colliding wind + power law model between 0.4 and lO.OkeV in 
the case of equal winds. The contributions of the three compo- 
nents are also individually displayed: the long-dashed, short- 
dashed, and dotted lines are respectively the mekal, wind- 
collision, and power law components. 

explore in the following section could be that the stars have 
winds of unequal strengths. 



7.2.2. Unequal winds 

Although Pachoulakis (1996 1 inferred only slightly different 
wind momenta (see Tabled, we note that the stellar parame- 
ters (T e ff , Lboi) derived by Pachoulakis yield very different val- 
ues of M (9 10~ 7 and 1 10~ 7 M yr _1 ) when passed through 
the relation of Vink et al. (2001 1. Such large differences in M, 
and hence wind strength, are however at odds with the essen- 
tially equal spectral types inferred from optical spectroscopy. 
Nevertheless, we attempted to establish whether an unequal 
wind CWB model could reconcile the predicted and observed 
luminosities. We thus developed a grid of models depending on 
four parameters: the terminal velocity of the primary (v Mi i), 
the mass loss rate of the primary (Mi), the terminal velocity 
of the secondary (Voo.2), and finally the wind momentum ratio 
defined as 

Ml Woo,! 

X] = - . 

M 2 Woo, 2 

This latter parameter is used to control the ratio of the wind 
strengths, to avoid any collision of the primary wind with the 
photosphere of the secondary. An additional normalization 
parameter as described for the equal wind case is also used. 
The velocity parameters are allowed to vary in a range between 
1700 and 2500 kms -1 , and the allowed values for Mi are 
taken between 3.0 10 -8 M Q yr -1 and 6.0 10 -6 M© yr -1 . The 
parameter r\ can vary between 1.0 and 1.44. Beyond this upper 
boundary the stronger wind of the primary crashes onto the 
secondary surface. 

The result of the fit seems only weakly sensitive to the value 
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of 77 (note that the range covered is quite narrow). If 77 is 
frozen to unity and velocities constrained to be equal, the same 
fitting result as described for the equal wind case is obtained. 
However, we were not able, with the current model, to obtain a 
satisfactory fit with a norm parameter close to unity. This result 
illustrates again the impossibility to reconcile the observed 
spectrum with the model at its current stage of development, at 
least within the parameter space we explored. 

8. Summary and conclusions 

The main results from our study of the XMM-Newton data of 
HD 159176 are the following: 

- The study of EPIC and RGS data reveals a soft spectrum 
which is consistently fitted with a two temperature thermal 
model. The hottest component is at about 0.6 keV for 
EPIC-pn and RGS, and about 1.0 keV for EPIC-MOS 
data. The X-ray luminosity is rather low, with a value of 
Lx/Lboi showing a moderate excess (~ 7) compared to 
what is expected for isolated stars with the same bolometric 
luminosity but without any wind-wind interaction. 

- Our analysis of RGS data reveals that lines are significantly 
broadened. For instance, the O VIII Ly a line at 19.0 A has 
a full width at half maximum of about 2500 kms -1 , in 
agreement with X-ray lines originating either from shocks 
distributed throughout the wind or from a colliding wind 
zone. 

- The study of EPIC light curves failed to reveal any 
significant variability on time scales of 100 - 2000 s. This 
indicates that, at least in a system like HD 159176, the 
hydrodynamic instabilities that might exist in the region 
of the shocked winds are not able to produce a clear 
variability of the X-ray emission. 

- The EPIC spectrum reveals a high energy tail which can 
not be fitted by thermal models (mekal or colliding 
wind). This hard X-ray emission component was fitted 
with a power law with a photon index of about 2.5. Its 
presence could reflect a non-thermal process such as 
inverse Compton scattering (Pollock ll987l Chen & White 
1991 ). The photon index is not too far from the 1.5 value 
expected for an X-ray spectrum arising from a population 
of relativistic electrons generated through an acceleration 
mechanism involving strong shocks. The possibility of 
a non-thermal X-ray component was already suggested 
for instance for a system like WR 110 by Skinner et al. 
( 2002 1. So far, the most prominent indication of relativistic 
electrons in stellar winds of early-type stars have been 
found in the radio domain where a significant fraction of 
the stars were found to display a non-thermal, probably 
synchrotron, emission (e.g. Bieging et al. 1989). In the case 
of HD 159176, radio observations failed to reveal such a 
non-thermal component, providing only an upper limit on 
the radio flux at 6 cm (Bieging et al. 1989). However, this 



does not rule out the possibility that relativistic electrons 
could be accelerated at the wind collision shock. In fact, 
the acceleration site would be buried so deeply within 
the radio photosphere that no synchrotron emission could 
escape and we would therefore observe HD 159176 as a 
thermal radio emitter. 

- Besides the non-thermal tail, the observed spectral shape 
can be consistently reproduced using a steady-state collid- 
ing wind model with a mass loss rate value of about 1.7 — 
2.6 10 -7 M Q yr _1 and a terminal velocity ranging between 
1850 and 1950 kms" 1 (or between 2140 and 2240 km s" 1 
for EPIC-MOS data). However, this model is unable to 
predict X-ray luminosities compatible with the observed 
spectrum of HD 159176. Theoretical values are systemat- 
ically higher than the observed X-ray luminosities. This 
disagreement between theory and observation is discussed 
hereafter. 

The disagreement between our XMM-Newton observation 
and the theoretical predictions could possibly be explained 
by several factors. First, let us recall that the kinetic power 
of the collision should be considered as an upper limit on 
the X-ray luminosity even in a highly radiative system. For 
instance, some of the collision energy might be taken away by 
the shocked gas which ends up with negligible velocity and 
pressure near the line of centers. In order not to pile up at the 
'stagnation point', this gas must be advected from the system, 
and the work needed to lift it out of the gravitational potential 
of the system could take away energy at the expense of X-ray 
emission. Second, higher values of the parameter 77 (see the 
unequal wind case, Sect. l7.2.2l should be considered. Although 
the current version of the model is unable to deal with the case 
where the primary wind crashes onto the secondary photo- 
sphere, this scenario should be envisaged. Let us emphasize 
however that the optical spectrum of HD 159176 does not pro- 
vide support for very large values of 77. Third, diffusive mixing 
between hot and cool material is likely to exist due to the 
instability of the shock front. As a consequence, the material 
tends to emit a much softer spectrum (i.e. EUV, or even UV), 
at the expense of X-rays. Unfortunately, current simulations 
lack the needed resolution to accurately deal with this mixing. 
Note that thermal conduction (Myasnikov & Zhekov 1998) 
is also expected to produce a softer spectrum. Finally, further 
developments of the current model are needed to address 
this issue. An improvement of the current steady-state model 
would be for instance to consider the effect of mechanisms 
able to lower the predicted luminosities like sudden radiative 
braking (Gayley et al. U997> . which can potentially occur 
every time unequal winds interact in a close binary system. In 
addition, orbital effects should also be included to study such 
systems. 

HD 159176 is the first short period colliding O + O binary 
studied with the high sensitivity of the instruments on board 
the XMM-Newton satellite. The major point of this system is 
that the shocks associated with the wind collision are radiative, 
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making them very difficult to simulate with current hydrody- 
namic models. For the first time, an alternative model has been 
used to address this case, following a steady-state geometri- 
cal approach leading to promising results. To achieve a better 
understanding of the special case of radiative colliding wind 
shocks, other close binary systems should be observed in con- 
junction with further developments of the new theoretical ap- 
proach followed in this study. 
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